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Excited B mesons Bi and are observed directly for the first time as two separate states in 
fully reconstructed decays to B'^^*\~ . The mass of Bi is measured to be 5720.6 ± 2.4 ± 1.4 MeV/c^ 
and the mass difference AM between and Bi is 26.2 ± 3.1 ± 0.9 MeV/c^, giving the mass of the 
B| as 5746.8 ± 2.4 ± 1.7 MeV/c^. The production rate for Bi and B| mesons is determined to be 
a fraction (13.9 ± 1.9 ± 3.2)% of the production rate of the B^ meson. 

PACS numbers: 12.15.Ff, 13.20.He, 14.40.Nd 



To date, the detailed spectroscopy of mesons contain- 
ing a b quark has not been fuUy established. Only the 
ground 0" states , B^, S", and the excited 1^ 
state B* are considered to be established by the PDG 
[H . Quark models predict the existence of two broad {Bq 
[JP = 0+] and BI [1+1) and two narrow {Bi [1+] and B* 
[2+]) bound P states [1 i, i, i, i, 0]. The broad states 
decay through an S wave and therefore have widths of a 
few hundred MeV/c^. Such states are difficult to distin- 
guish, in effective mass spectra, from the combinatorial 
background. The narrow states decay through a D wave 
and therefore should have widths of around 10 MeV/c^ 
[3i @, H- Almost all observations of Bi and B2 have 
been made indirectly in inclusive or semi-inclusive de- 
cays d, H, U I , which prevents their separation and 
a precise measurement of their properties. The measure- 
ment by the ALEPH collaboration fl2\, although par- 
tially done with exclusive B decays, was statistically lim- 
ited and model dependent. The masses, widths, and de- 
cay branching fractions of these states, in contrast, are 
predicted with good precision by various theoretical mod- 
els 0, S 0, H, IE 01 ■ These predictions can be verified 
experimentally, and such a comparison provides impor- 
tant information on the quark interaction inside bound 
states, aiding further development of non-perturbative 
QCD. This Letter presents a study of narrow L = 1 states 
decaying to i3+(*'7r~ with exclusively reconstructed 5+ 
mesons using data collected by the DO experiment during 
2002-2006 and corresponding to an integrated luminos- 
ity of about 1.3 fb^^. Throughout this Letter, charge 
conjugated states are implied. 

The DO detector is described in detail elsewhere [l3| . 
The detector components most important for this anal- 
ysis are the central tracking and muon systems. The DO 
central tracking system consists of a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both 
located within a 2 T superconducting solenoidal mag- 
net, with designs optimized for tracking and vertexing 
at pseudorapidities \?]\ < 3 and |ry| < 2.5, respectively 
(where 77 = — ln[tan(6'/2)], and 9 is the polar angle mea- 
sured with respect to the beam line). The muon system 
is located outside the calorimeters and has pseudorapid- 
ity coverage jryj < 2. It consists of a layer of tracking 
detectors and scintillation trigger counters in front of 
1.8 T toroids, followed by two similar layers behind the 
toroids 



decays: 



The Bi and B* 
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B+TT- 
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mesons are studied by examining 



J5"'"(*V candidates. This sample includes the following 



The direct decay Bi — > B'^tt~ is forbidden by conserva- 
tion of parity and angular momentum. The B'^ meson 
is reconstructed in the exclusive decay B^ J/tpK^ 
with J/ip decaying to fi'^ ii~ . Each muon is required to 
be identified by the muon system, have an associated 
track in the central tracking system with at least two 
measurements in the SMT, and a transverse momentum 
pt^ > 1.5 GeV/c. At least one of the two muons is re- 
quired to have matching track segments both inside and 
outside the toroidal magnet. The two muons must form 
a common vertex and have an invariant mass between 
2.80 and 3.35 GeV/c^, to form a J/ip candidate. An 
additional charged track with px > 0.5 GeV/c, with to- 
tal momentum above 0.7 GeV/c and with at least two 
measurements in the SMT, is selected. This particle is 
assigned the kaon mass and required to have a common 
vertex, with < 16 for 3 degrees of freedom, with the 
two muons. The displacement of this vertex from the pri- 
mary interaction point is required to exceed three stan- 
dard deviations in the plane perpendicular to the beam 
direction. The primary vertex of the pp interaction was 
determined for each event using the method described 
in Ref. 15]. The average position of the beam-collision 
point was included as a constraint. 

From each set of three particles fulfilling these require- 
ments, a candidate is constructed. The momenta of 
the muons are corrected using the J/tp mass [H as a con- 
straint. To further improve the B~^ selection, a likelihood 
ratio method [l6| is utilized. This method provides a sim- 
ple way of combining several discriminating variables into 
a single variable with increased power to separate signal 
and background. The variables chosen for this analy- 
sis include the smaller of the transverse momenta of the 
two muons, the of the i?+ decay vertex, the B~^ de- 
cay length divided by its error, the significance (defined 
below) Sb of the B'^ track impact parameter, the trans- 
verse momentum of the kaon, and the significance Sk of 
the kaon track impact parameter. 

For any track i, the significance Si is defined as Si — 
y^[eT/cr{eT)]'^ + [eL/(T{eLW, where er (el) is the projec- 
tion of the track impact parameter on the plane perpen- 
dicular to the beam direction (along the beam direction), 
and a{eT) is its uncertainty. The track of each 
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FIG. 1: Invariant mass distribution of J/il^K events. The 
solid line shows the sum of signal and background contribu- 
tions, as described in the text. The contribution of J/ip-K^ 
events is shown by the solid filled area and the J /ipK^* contri- 
bution is shown by the hatched area. The dashed line shows 
the exponential function modeling the combinatorial back- 
ground. 



B'^ is formed assuming that it passes through the recon- 
structed vertex and is directed along the reconstructed 
momentum. 

The resulting invariant mass distribution of the 
J/^K~^ system is shown in Fig. [TJ The curve represents 
the result of an unbinned likelihood fit to the sum of 
contributions from B~^ J/tpK^, i?+ J/tpir^, and 
5+ — !■ J/tpK^* decays, as well as combinatorial back- 
ground. The mass distribution of the J/ijjK^ system 
from the B~^ J/^K~^ hypothesis is parameterized by 
a Gaussian with the width depending on the momentum 
of the K'^. For the contribution from B^ J/ijjTT^ 
decays, the width of the J/^tt^ mass distribution is 
parametrized with the same momentum-dependent width 
as the B+ J /ifjK^ decays, and then transformed to 
the J/tpK^ system by assigning the kaon mass to the 
charged pion. The decay B J/ij^K^* with Kn 
produces a broad J/ipK^ mass distribution with the 
threshold near M{B) — M(7r). It is parameterized using 
Monte Carlo simulation (described later). The combina- 
torial background is described by an exponential func- 
tion. The 5+ — > J/'0i^+ and i?+ — > J/'07r+ mass peaks 
contain 23287 ±344 (stat.) events. 

For each reconstructed B meson candidate with mass 
5.19 < M{B+) < 5.36 GeV/c^, an additional charged 
track with transverse momentum above 0.75 GeV/c and 
charge opposite to that of the B meson is selected. The 
selection 5.19 < M{B+) < 5.36 GeV/c^ reduces the 
number of B~^ candidates to 20915 ± 293 (stat.). Since 
the Bj mesons (where Bj denotes both Bi and i3|) de- 
cay at the production point, the additional track is re- 
quired to originate from the primary vertex by applying 
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FIG. 2: Invariant mass difference AM = M{B+-k-) - M{B+) 
for exclusive B decays. The line shows the fit described in the 
text. The contribution of background and the three signal 
peaks are shown separately [color online]. 



the condition on its significance < V6- 

For each combination satisfying the above criteria, the 
mass difference AM = M(B+7r") - M{B+) is com- 
puted. The resulting distribution of AM is shown in 
Fig. [21 The signal exhibits a structure that is inter- 
preted in terms of the decays ([TH3]). Since the pho- 
ton from the decay B* — > Bj is not reconstructed, 
the three decays should produce three peaks with cen- 
tral positions Ai ~ M{Bi) — M(B*), corresponding to 
the decay Bi B*n, A2 = M{B^) - M{B*), corre- 
sponding to B^ B*TT, and A3 = M{B^) - M{B), 
corresponding to B2 Btt. Note that in this case, 
A2 = A3 - [M{B*+) - M{B)] = A3 - 45.78 MeV/c^ 
[l|. Following this expected pattern, the experimental 
distribution is fitted to the following function: 

F{AM) = ^^,ig(AM)-HFbckg(AAf), 
i^sig(AM) = N-{fi-D{AM,Ai,ri) 

+ {l-fi)-[f2-D{AM,A2,T2) 

+ (l-/2)-7^(AAf,A3,r2)]}. (4) 

In these equations, Fi and r2 are the widths of Bi and 
B2, fi is the fraction of Bi contained in the Bj signal, 
and /2 is the fraction of B2 B*tt decays in the B2 sig- 
nal. The parameter N gives the total number of observed 
Bj _B+(*V decays. The background Fbckg(AAf) is pa- 
rameterized by a fourth-order polynomial. 

The function D{x, xo,T) in Eq. (|3|) is the convolution of 
a relativistic Breit-Wigner function with the experimen- 
tal Gaussian resolution in AM. The width of resonances 
in the Breit-Wigner function takes into account threshold 
effects using the standard expression: 



T{x) = Fo— — 



(5) 
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The variables k and fcg in Eq. ([5]) are the magnitudes 
of the pion three-momentum in the Bj rest frame when 
Bj has a four-momcntum-squared equal to x'^ and Xq, 
respectively. Tq is the total decay width, and F'-'^^kj kn) 
is the Blatt-Weiskopf form factor for L = 2 decay [l|, ll7| . 

The resolution in AM is determined from simulation. 
All processes involving B mesons are simulated using 
the EvtGen generator [l8| interfaced with pythia [l^, 
followed by full modeling of the detector response with 
GEANT [i^] and event reconstruction as in data. The dif- 
ference between the reconstructed and generated values 
of AM is parameterized by a double Gaussian function 
with the a of the narrow Gaussian set to 7.5 MeV/c^, the 
a of the wide Gaussian set to 17.6 MeV/c^, and the nor- 
malisation of the narrow Gaussian set to 3.8 times that 
of the wide Gaussian. Studies of the — > J/ipK^ mass 
peak show that simulation underestimates the mass reso- 
lution in data by ~10%. As such, the widths of the Gaus- 
sians which parameterise the Bj resolution are increased 
by 10% to match the data, and a 100% systematic uncer- 
tainty is assigned to this correction. The widths of the 
observed structures are compatible with the experimen- 
tal mass resolution, and the fit is found to be insensitive 
to values of Fi and r2 below the mass resolution with 
the current statistics. Therefore, both these widths are 
fixed at 10 MeV/c^ in the fit, as suggested by theoretical 
models 0, 0] ■ They are varied together over a wide range 
to estimate the associated systematic uncertainty. 

With these assumptions, the following parameters of 
Bi and i?2 are obtained: 



M{Bi) - M{B+) 
M{Bl) ~ M{Bi) 



441.5 ±2.4 ±1.3 MeV/c^ 
26.2 ±3.1 ±0.9 MeV/c^. (6) 



where the first uncertainty is statistical, and the second 
is systematic. The correlation coefficient of these mass 
measurements is —0.659. With these relations, and using 
the mass oftheB+ (H, the absolute masses of the Bi and 



B* 



are: 



M{Bi) = 5720.6 ± 2.4 ± 1.4 MeV/c^ 
M{B;) = 5746.8 ± 2.4 ± 1.7 MeV/c^ 



(7) 



The number of Bj decays is found to be iV = 662 ± 91. 
The xVd.o.f. of the fit is 33/40. Without the Bj signal 
contribution, the x^/d.o.f. of the fit increases to 97/45, 
which implies that this structure is observed with a sta- 
tistical significance of more than 7a. Fitting with only 
one peak increases the x^/d.o.f. to 54/42, which corre- 
sponds to more than a 4(T significance that more than 
one resonance is observed. With the B^ — *■ B*^Tr decay 
removed from the fit, the x^/d-o.f. of the fit increases to 
41/41. Although with the current statistics we can not 
distinguish between the two- and three- peaks hypothe- 
ses, theory suggests that i3| decays with almost equal 
branching ratios into Btt and B*~^Tr [3, Q, and our fit 
indeed indicates a preference for this expected pattern. 



The number of Bj mesons and values /i and /2 ob- 
tained from the fit are used to measure the production 
and decay ratios of Bi and B2: 

BrjBi ^ B*+7t) eo 
■Ki - „ T^Tzr^ - h - — 



i?2 = 



Br{Bj B(*)tt) 
= 0.477 ± 0.069 ±0.062 
Br{B^ B*tt) 



= f2-- 



R.J = 



Br{B* 
= 0.475 ± 0.095 ±0.069, 

Br{b B°^ B(*)7r) 



3-NiBj 



Br{b B+) 
= 0.139 ±0.019 ± 0.032. 



2-NiB+)-eo 



(8) 



Here ei, £2 and £3 are the efficiencies to select an 
additional pion from the Bj decay for decay modes 
Bi B*+TT-, BI B*+TT- and B^ B+tt", re- 
spectively. They are determined from a simulation sepa- 
rately for each decay mode ([l] [3]). The overall efficiency 
for detecting a pion from any Bj — > _B+(*'7r~ decay is 
£0 = 0.342 ± 0.008 ± 0.028. The value for Rj takes into 
account the decay Bj _B°7r° assuming isospin conser- 
vation. 

For the Bj mass fit, the infiuences of different sources 
of systematic uncertainty are estimated by examining the 
changes in the fit parameters under a number of varia- 
tions. Different background parameterizations are used 
in the fit to the AM distribution. In addition, the ef- 
fect of binning is tested by varying the bin width and 
position. The parameters describing the background are 
allowed to vary in the fit and their uncertainties are in- 
cluded in our results. To check the effect of fixing Fi and 
F2 at 10 MeV/c^, a range of widths from to 20 MeV/c^ 
is used. The effect of the uncertainty on the mass dif- 
ference M{B+*) - M{B+) [H is also taken into account. 
Different parameterizations of the detector mass resolu- 
tion are tested, and in addition the fit is made without 
the 10% mass resolution correction. The uncertainty in 
the absolute momentum scale, which results in a small 
shift of all measured masses, is also taken into account. 
The summary of all systematic uncertainties in the Bj 
mass fit is given in Table [H 

The measurement of the relative production rate Rj 
uses the pion detection efficiencies predicted in simula- 
tion, as well as the numbers of Bj and B~^ events. To 
estimate the systematic uncertainty on the number of 
5+ events, different parameterizations of the signal and 
background are used for the fit. The resulting uncer- 
tainty is ±200 B~^ events. The systematic uncertainty 
on the number of Bj events is ±140 (see Table The 
uncertainty of the impact parameter resolution in the 
simulation is estimated to be «10% jl^]. It can infiuence 
the measurement of the selection efficiency of the pion 
from the Bj decay, and its contribution to the system- 
atic uncertainty of Rj is found to be 0.0056. The track 
reconstruction efficiency for particles with low transverse 
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TABLE I: Systematic uncertainties of the Bj parameters determined from the AM fit. The rows show the various sources 
of systematic error as described in the text. The columns show the resulting uncertainties for each of the five free signal 
parameters as described in Eq. Q. AM{Bi) and A[M(B|) - M{Bi)] are in MeV/c^. 



Source 


AA/(i3i) 




A_Ri 


AR2 


AA'' 


Background parameterization 


0.15 


0.15 


0.010 


0.009 


19 


Bin widths/positions 


0.85 


0.70 


0.006 


0.026 


12 


Value of r 


0.75 


0.55 


0.023 


0.032 


138 


B"*"* mass uncertainty 


0.30 


0.25 


0.004 


0.004 


6 


Momentum scale 


0.50 


0.03 


0.000 


0.000 





Resolution uncertainty 


0.20 


0.05 


0.007 


0.004 


10 


Efficiency uncertainties 






0.056 


0.054 




Total 


1.30 


0.90 


0.062 


0.069 


140 



momentum is measured in Ref. 21;] and good agreement 
between data and simulation is found. This comparison 
is valid within the uncertainties of branching fractions of 
different B semileptonic decays, which is about 7%. This 
uncertainty results in a 0.0096 variation of Rj. An ad- 
ditional systematic uncertainty of 0.0008 associated with 
the difference in the momentum distributions of selected 
particles in data and in simulation is taken into account. 
Combining all these effects in quadrature, the total sys- 
tematic uncertainty in the relative production rate Rj is 
found to be 0.032, of which the dominant contribution 
comes from the uncertainty on the number of Bj events. 

Different consistency checks of the observed signal are 
performed. The stability of the fit under different se- 
lections of TT meson from Bj decay and the 
width of the B^ mass window is verified. Events with 
positively and negatively charged pions are analyzed sep- 
arately, and consistent results are obtained. A comple- 
mentary sample of events containing a pion not compati- 
ble with the primary vertex is selected, and no significant 
Bj signal is observed. Events with wrong charge combi- 
nations (_B+7r+ and B^tt^) also show a signal consistent 
with zero. In addition, the fit is repeated without the 
Blatt-Weiskopf form-factor and no visible change in re- 
sults is observed. 

In conclusion, the Bi and i?2 mesons are observed for 
the first time as two separate states. Their measured 
masses are given by Eq. ([6]). The Bj production rate, 
the branching fraction of B2 to the excited state B*, 
and the fraction of the Bi meson in the Bj production 
rate are also measured as given in Eq. ([8]). These results 
will help to develop models describing bound states with 
heavy quarks. 
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